Below you can find a link to the MATLAB code that calculates the L2 bitmap projection using the isogeometric finite element method. In particular, this code implements a variable-direction solver algorithm.

The code execution is also possible in the free Octave environment.
Download code or see Appendix 4.

We set the variables with the path to the input file in tif format
\( filename = 'C://Users/Maciej/Dropbox/bitmapa.tif' \)
then we give the number of mesh elements in the x and y directions, and the degrees of the B-spline function in these directions
\( nx=4 \)
\( ny=4 \)
\( px=2 \)
\( py=2 \)
Then we start the first procedure
\( bitmap\_splitting(filename,nx,ny,px,py) \).
The code, after building the system of equations and solving it, draws the bitmap projections in the open window.

Exercise 1: Projection of the terrain

Content of exercise:

Please convert the bitmap projection code to draw projections of the area.

Solution:

Solution: download code or see Appendix 4A.